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1 Introduction 



Many efforts in contemporary pliysics are devoted to tlie study of quantum 
systems at finite temperature or density, like e.g. particles propagating in hot, 
compressed nuclear matter. While the vacuum state of quantum field theory 
is invariant under Lorentz transformations, matter states have, in general, 
less symmetry. This leads to a modification of the particle spectrum needed 
for the description of such a state 0. For the application of perturbative 
(i.e. diagrammatic) methods of quantum field theory, this modification has 
two main effects. 

First, one has to use 2x2 matrix valued propagators, as in Thermo 
Field Dynamics (TFD)0 or the Closed-Time-Path formalism [Q. Second, 
a consistent description requires elementary excitations with a continuous 
mass spectrum, rather than physical quasi-particles of infinite lifetime 

A well defined diagrammatic expansion at finite temperature therefore 
takes into account a spectral function of "particles" which deviates from a 
^-function, i.e., which is more than a mass-shell constraint. In several papers 
it was shown recently, how the propagators of such an expansion can be 
diagonalized, and what the physical meaning of the matrix structure is in 
the framework of Thermo Field Dynamics ^ • 

In the present work, a calculational recipe for TFD is discussed, which 
was developed for non-relativistic theories ^ . The method is extended to rel- 
ativistic quantum fields, where it makes the calculation of Feynman diagrams 
for matrix valued propagators a very simple task. To this end, a quantum 
field theory with minimal coupling of fermions and scalar uncharged bosons 



in thermal equilibrium is considered (see refs. |TOl for applications). 

In TFD, the thermal instability of observable states is absorbed into a 
Bogoliubov transformation. This Bogliubov transformation can be written 
also for interacting systems, where it defines stable, albeit non-observable 
quasi-particles |^, |^, |^. At the level of the quasi-particle operators, it can 
be cast into a single matrix form for the bosonic and the fermionic sector of 
the model, with a parameter n that resembles the phase-space occupation 
factor. While one has some freedom in parameterizing the transformation, 
the most useful way to do so is 



BsAn) = , (1) 




2 



since it is linear in n [Q]. In the fermionic sector, the n are Fermi-Dirac 
functions, labeled rip for positive energy states and rip for negative energy 
states. The thermal quasi-particle picture of TFD requires that they depend 
on a continuous energy variable. 

In the bosonic sector, n is a Bose-Einstein distribution function. 



2 Full propagators in thermal equilibrium 

The full bosonic two-point Green's function, or propagator, at finite temper- 
ature has been derived in refs. 0, as 

oo 

D^^''\ko,k) = jdEpB{E,k) X (4) 



ko—E—ie 
1 \ 



- r,Bl{nB{E))[ >'o^^- ^_ {Bl{nB{E)))-' . 

\ ko+E+ie / J 

Pb{E, k) is a positive weight function with support only for positive energy 
arguments, = diag(l, —1). 

In the above equation, the terms propagating particle and anti-particle 
states have been kept separately. Since eqn. can be continued to negative 
energy arguments, the above propagator simplifies to 

oo 

DH)(A;o,fc) = JdEABiE,k) X (5) 



{BBiriBiE)))-' [ ^«-^+" 1 BBiuBmrs. 

\ ko-E-ie / 
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AsiE, k) is the spectral function of the boson field, 

Ab{E, k) = pb{E, k)Q{E) - pb{-E, k)e{-E) , (6) 

and the limit of free particles with mass m is recovered when 

AB{E,k) — ^ sign(^) 5(^2 - - m^) . (7) 

The fermion propagator is calculated in a similar way. The occurence of a 
nonzero chemical potential does not hinder the above simplification, since 
the occupation number parameters for positive and negative energy states 
are related by 

1 - n^{-E) = nUE) ^ nj.{E) , (/^ (-i?))"^ = /^(E) ^ ME) . (8) 

Some elementary matrix algebra then gives, in terms of the fermionic spectral 
function Af{E,p) which has support on the whole real energy axis, 

oo 

S'k'\po,p)= JdEAp{E,p) X 

— oo 

niBp{np{E)))-' (^^^^ 1 ) BpME)). (9) 

\ po—E—ie / 

The limit of free Dirac particles with mass M is recovered when 

Af{E, p) {p^Y + M) sign(^) 6iE^ - p" - M^) . (10) 

Note, that the above equations for the full propagators contain the full con- 
tribution of the negative energy states. 



3 Column vector notation for propagators 

To simplify the use of such matrix valued propagators, a column vector no- 
tation was introduced in ref. 0. It is obtained by rewriting the thermal 
Bogoliubov matrices defined in eqn. (|l|) as the outer product of two column 
vectors according to 
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This decomposition works only for matrices with zero determinant. The art 
in applying this idea therefore lies in splitting the propagator matrices into 
a sum of terms, each with determinant zero. 

For the above propagators this task is trivial: one simply isolates the 
parts propagating forward and backward in time, this gives the full fermion 
and boson propagators as 

(12) 

Po-E-ie{f^\E)}{l}) ' 



—oo 





D(«'')(A;o,fc) = JdEAB{E,k)nB{E)x (13) 

1 [ 1 

k,-E-ze\fB\E) 

Note, that each column vector is associated with one endpoint of a propagator 
line. Consequently, an A^-point function of the interacting theory is described 
as an object composed of N column vectors, 

r'"' = (?U?|--(^l . (14) 



Now consider a typical diagram of the perturbative expansion of an inter- 
acting quantum field theory. In such a diagram, in general bare vertices and 
irreducible vertex functions are joined by propagators. The rule for combi- 
nation of their column vectors now is the following: 

Rule 1. If two or more propagator lines, counted by an index 1 < I < L, 
meet in one L-point vertex, and the vertex is still "open", i.e., external lines 
can be attached to it, combine their column vectors to a new column vector 
as 

{::}^{:}--{n;;}- 
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Rule 2. If the vertex is "closed", i.e., if it is an interior vertex and the last 
possible propagator line with index L is attached to it, contract it to a scalar 
function 




Note, that in ref. |^ a relative negative sign was used in this rule. In the 
present work, this was absorbed into the definition of the coupling constant, 
see below. 



4 Applications of the method 

In a first step, the above rules are applied to "one-loop" self energy diagrams. 
Since the propagators of the model contain the full spectral functions, this 
"one-loop" approximation can go substantially beyond a naive perturbation 
theory. First, consider the one- loop polarization tensor, 

where a, b are the thermal matrix indices taking values 1 and 2. The coupling 
constants are (cf. ref. 0) (7^ = = g, and the exchange of the indices in 
the second factor is achieved by a simple transposition of the corresponding 
column vectors. Hence one has 

,4 00 

n('^'')(A;o,fe) = -^ JdEdE'Ti[AF{E,p + k)AFiE',p)] nF{E)nF{E')x 

— 00 

[po + ko-E + ze{l}{^''-l Po + ko-E-ie{fF\E)}{l}) 

[po-E' + zei^'^-l p,-E'-ze{l}{fF\E')}) 

where the last x indicates that the product is closed, i.e. the trace also 
applies to the thermal indices. 

In performing the po-integration, out of the four possible products in 
the integrand only those survive which have different relative sign of ie in 



X 



X 
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their factors. Applying rule 1, the complete polarization tensor therefore is 
obtained as 

— CO 

( 1 [-!][ f-AE) \ 1 \rF\E')\ \-l\\ 

\k, + E'-E + ie\ 1 /\-//(i?')J h + E'-E-te\-fp\E)]\ 1 ]) ' 

In the next step one uses, that this self-energy function appears as the kernel 
of the Schwinger- Dyson equation for the full boson propagator as 

D{k) = Do{k) + DQ{k)»Y{{k)»D{k) (20) 

To pick the retarded full propagator out of this sum, keep only the terms 
with positive ie in the denominator (the retarded propagator is an analytical 
function in the upper complex fco-plane). Thus, the second term on the r.h.s 
of this equation has, in the retarded part, the column vector structure 

In other words, inserting an arbitrary number of polarization functions re- 
tains the same column vector structure as in eqn. (|T^) ! One can thus read 
off the above equations the expression for the retarded polarization function, 
which completely determines the retarded full boson propagator, as 

— oo 

(22) 

In ref. [jlO[, this equation was made the basis of an explicit calculation of a 
bosonic spectral function at finite temperature. 
The second example is the Fock diagram, 

S('^^)(po,p)=^/^Tr[(?('^)5('^^)(po + fco,P + %('^/^('"nfco,fc)] . (23) 
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Following the same steps as above, i.e., contracting the column vector prod- 
ucts at both vertices, and performing the /co-integration such that only terms 
with different sign of ie survive, then gives 

,3, 00 

S^'^'HpcP) = Jj^ I dEdE' AFiE,p + k)ABiE',k) nF{E)nB{E')x 

— 00 

[p, + E'-E + te { -1 } { 11\e\ ] + ^;VW^¥^e [ ^{e) } { -1 }) (^4) 

Again, put this into a Schwinger-Dyson equation 

S{p) = So{p) + So{p)*np)*S{p), (25) 

and isolate the retarded full fermion propagator by keeping only the terms 
with positive ie in the denominator. The second term on the r.h.s then has, 
in the retarded part, the factors 

nAPo) 1 1 } I ^^l^f ^ } {riBiE') - . (26) 

Again, the column vector structure of the full propagator (|T2|) is conserved by 
the insertion of self energy diagrams. From the above expression furthermore 
follows the retarded self energy function as 

S^(P0,P)=/ JdEdE'AAE,p+k)As{E',k) Mfl'l^^f^ ] • 

—00 

(27) 



5 Conclusion 

In the present work, the column vector calculus of Thermo Field Dynamics 
was introduced for relativistic quantum fields. The simplicity of its applica- 
tion was shown in two examples: two well known equations for retarded self 
energy functions have been obtained within a few lines. 
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While this might not be a great gain for the generahzed one-loop examples 
given here, it makes higher order calculations much simpler than with other 
methods. An example for a non-relativistic vertex correction is dicussed in 
ref. p[. The conclusion is, that the column vector calculus is even simpler 
to apply than the Matsubara technique. 



In ref. |]T0| , the above eqn. ( P2|) was used for an application of TFD. There 
we considered a model with two different types of fermions, i.e., nucleons and 
A33 resonances, coupled to a pion field, and calculated the pionic spectral 
function. More sophisticated applications are in progress. 
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